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^ ■ Abstract 
<N 

The first differential cross section for Mott scattering of a Dirac-Volkov electron is 
reviewed. The expression (26) derived by Szymanowski et al. [Physical Review A 56, 
O . 3846,(1997)] is corrected. In particular, we disagree with the expression of (^) they 

obtained and we give the exact coefficients multiplying the various Bessel functions 

■ appearing in the scattering differential cross section. 

o ' 
O 

a: 

p_; 1 Introduction 

In a pioneering paper, Szymanowski et al. [Q] have studied the Mott scattering process 
in a strong laser field. The main purpose was to show that the modifications of the Mott 
scattering differential cross section for the scattering of an electron by the Coulomb potential 
of a nucleus in the presence of a strong laser field, can yield interesting physical insights 
concerning the importance and the signatures of the relativistic effects. Their spin dependent 
relativistic description of Mott scattering permits to distinguish between kinematics and spin- 
orbit coupling effects. They have compared the results of a calculation of the first Born 
differential cross section for the Coulomb scattering of the Dirac-Volkov electrons dressed by 
a circularly polarized laser field to the first Born cross section for the Coulomb scattering of 
spinless Klein-Gordon particles and also to the non relativistic Schrodinger-Volkov treatment. 
The aim of this comment is to provide the correct expression for the first-Born differential cross 
sections corresponding to the Coulomb scattering of the Dirac-Volkov electrons. One the one 
hand, We show that the terms proportional to sm(20o) are missing in [1J], where 0o is the phase 
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stemming from the expression of the circularly polarized electromagnetic field. The claim of 
I] that they vanish is not true. These terms do not depend on the chosen description of the 
circular polarization in cartesian components. On the other hand, We perform the calculations 
with some details and throughout this work, we use atomic units (h = e = m = 1) where m 
denotes the electron mass. The abbreviation DCS stands for the differential cross section. 

The organization of this paper is as follows : in Section 2, we establish the expression of the 
S'-matrix transition amplitude as well as the formal expression of scattering DCS. In Section 
3, we give a detailed account on the various trace calculations and show that indeed there is 
a missing term proportional to sin(20o) that is not equal to zero. This term as well as a term 
proportional to cos(20 o ) contribute to (^) and multiply the product J s+1 (z) J s ^i(z) where 
J s {z) is an ordinary Bessel function of argument z and index s. The argument z appearing 
in the above mentioned product will be defined later. Then, we carry out the derivation of 
the correct expression of the scattering DCS associated to the exchange of a given number of 
laser photons. We end by a brief a conclusion in Section 4. 



2 The S'-matrix element and the scattering differential 
cross section. 

Exact solutions of relativistic wave equations are very difficult to obtain. However, in 
seminal paper, Volkov || obtained the formal solution of The Dirac equation for the relativis- 
tic electron with 4-momentum p M inside a classical monochromatic electromagnetic field A^. 
These solutions are called the relativistic Volkov states. The plane wave electromagnetic field 
of 4-momentum k^ {k^ = k 2 = 0) depends only on the argument M and 

therefore A^ is such that : 

A 1 * = A»(k.x) = A"(</>) (1) 

The 4-vector A^ satisfies the Lorentz gauge condition d^A^ = or equivalently k^A^ = 0. 
The Dirac- Volkov equation in an external field A^ is: 

| (p - - c Af - c 2 - y/^ U } ^) = ( 2 ) 

where F^ u is the electromagnetic field tensor F^ y = d^A u — d v A^ and o^ v = |[7 M ,7^]. The 
matrices 7^ are the anticommuting Dirac matrices such that 7^7^ + 7^7^ = ig^X^ where g^ u 
is the metric tensor g^ v = diag(l, —1,-1,-1) and I4 is the identity matrix in four dimensions. 
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The solutions of Eq. (||]) are the relativistic Dirac-Volkov wave functions: 

il>Jx) = R(p)?$=l£Le lS ^ (3) 

where : 



and the function S(x) is given by : 

S( x ) = -p.x - / — — 

Jo C l ft -J 



d£ (5) 



c(k.p) 

In Eq. (^|), u(p,s) represents a Dirac bispinor which satisfies the free Dirac equation and 
is normalized according to u(p, s)u(p, s) = u*(p, s)^°u(p, s) = 2c 2 . We consider a circularly- 
polarized field : 

A = ai cos(0) + Gt2 sin(0) (6) 

where <fi = k.x. We choose a\ = a\ = a 2 = A 2 and a\.a2 = g^-^i = 0. The Lorentz condition 
k.A = implies a\.k = a^.k = 0. If one assumes that A^ is quasi-periodic so that its time 
average is zero A^ = 0, then using the Gordon identity, the averaged 4-current is easily 
obtained : 

If one sets : 



this yields : 



with : 



p 1 2c 2 (k.p)' 

Q-Q = = ™\<? (9) 

ml = l-- (10) 

One often calls the averaged 4- momentum g M a quasi-impulsion. Note that g M = (Q/c,q). 
The quantity m» plays the role of an effective mass of the electron inside the electromagnetic 

field. For the study of the process of Mott scattering in presence of a laser field, we use the 
Dirac-Volkov wave functions @ normalized in the volume V: 



where 



Rip) = R(q) = i + U = 1 + ^(kpj^ 1 cos(0) + ^ 2 sin(0)) (12) 
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and : 

5(9, x) = -^-—an^—cosM 

- ---M Sm(0) + M C ° S(0) (13) 
We turn now to the calculation of the transition amplitude. The interaction of the dressed 
electrons with the central Coulomb field : 

^ = (A,0) (14) 

is considered as a first-order perturbation. This is well justified if Za <C 1, where Z is the 
nuclear charge of the nucleus considered and a is the fine-structure constant. We evaluate the 
transition matrix element for the transition (i — > /) : 

S fi = — [ d 4 xip qf (x)j—ip qi (x) (15) 
C J ' \x\ 

We first consider the quantity : 

= T^^TIp^' *f)RM^R(Pi)»<Pi, 8i )e-< 8 ^- 8 ^ (16) 

We have : 

e -i(S( qf ,x)-s( qi ,x)) = eX p[i( g/ - qi ). x - izsm((f) - O )] (17) 
where z is such that : 

z = jal + al (18) 
whereas the quantities a.\ and a 2 are given by: 

^ = (aj.pj) _ (ai.pf) ^ = (a 2 .pi) _ (a 2 -P/) ^ 
c(k.pi) c{k.pfY c(k.pi) c(k.pf) 

and the phase 0o is such that 0o — arccos(a:i/,2) = arcsin^/z) = arctan(o;2/«i)- It is 
important at this stage to perform intermediate calculations in order to reduce the numbers 
of 7 matrices that will appear when one calculates the scattering DCS. After some algebraic 
manipulations, one gets : 

u(p f , Sf)R(p f )^R{pi)u{p h = u(p f , s f ) [C + Ci cos(0) + C 2 sin(0)]w(pi, s,) (20) 

where the three coefficients Co, C\ and C 2 are respectively given by : 

C = 7° - 2k a 2 MPiXPf) 

C 2 = c( Pi )^U 2 + c(p f )4 2 ^° (21) 
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with c(p) = 2c ^ k p j and ko = k° = uj/c. Therefore, the transition matrix element becomes : 



Sfi = % — J d 4 x -j=^=-j== u{p h s f )[C + d cos(0) + C 2 wn((/>)]u(pi, Si) 

x exp[i(qf — qij.x — izsin{(f> — (f> )] (22) 

We now invoke the well-known identities involving ordinary Bessel functions J s (z) : 

i ) oo r b s ) 

cos(0) I e -«™(*-*>) = ^ J ^ I e -i^ (23) 
sin(0) I s=-oo I £ 2s I 



with : 









hi 


I bL J 





J s (z)e is<t>0 

(J s+1 (z)e^ s+1 ^° + J s _ 1 (z)e i ^ s - 1 ^ )/2 J> (24) 
(J s+ i(z)e i(s+1) *° - J s - 1 (z)e i ( s - 1 ^°)/2i 



Evaluating the integrals over x and x yields for Sfi : 



= i47rZ y> 27r£(Q / -Q i -sa;) (a) 

/l v^Fv^V s f4 \q f - qi -8k\* fi 1 J 

where the quantity Mj/ is defined by : 

Ki = "(P/> s/MCbS, + Ci^i- + C^Mp*, Si ) (26) 

To evaluate the DCS, we first evaluate the transition probability per particle into final states 
within the range of momentum dq^ : 

dW ft - \S fi \ j-^ 

(4vr) 2 Z 2 ^ T 2-n5{Q f -Q l -sw) (s)2 Vdq f 
2Ql V.2Q f V \ qf - q .-sk\* 1 ^'(2^ 1 ] 

where we have used the rule of replacement : 

[2nS(Q f - Qi - sw)] 2 -> 2n5(0)2n5(Q f - Q t - sw) 

= T2n5(Q f -Qi- sw) (28) 

Next, we have for the transition probability per unit time : 

dR " ~ — ~ 2 Qi V2 Ql v 2- | 9/ |M /- 1 (So? (29) 
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Dividing dRji by the flux of incoming particles : 



\Qi\c z 



(30) 



then using the relation |<fy|cf|<fy| = jsQfdQf and integrating over the final energy, we get for 
the scattering DCS : 



da Z 2 \q } \ 



dn f c 4 \ Qi 



E 



|M 



(»)|2 



fi 



q f -Qi- sk\ 



E 



da^ 



Qf=Qi+sw 



dtt. 



(31) 



Qf=Qi+sw 



where : 



da^ 



dtt 



f 



Z*\q f \ 



\M 



(*)|2 



fi 



\Qi\ \Qf-Qi~sk\ 



(32) 



Q f=Qi+sw 



Qf=Qi+sw 

The calculation is now reduced to the computation of traces of 7 matrices. This is routinely 
done using Reduce ||. We consider the unpolarized DCS. Therefore, the various polarization 
states have the same probability and the actually measured DCS is given by summing over the 
final polarization Sf and averaging over the initial polarization Sj. Therefore, the unpolarized 
DCS is formally given by : 

da A da^ 



dtt 



f 



E 



dtt 



f 



Qf=Qi+sw 



where 



da^ 



dtt, 



Qf=Qi 



z 2 


9/1 


1 


c 4 


Qi\ 


\Qf -Qi-sk 


1 



EEKi 2 



(33) 



(34) 



Qf=Qi+sw 



3 Trace calculations. 



Since the controversy is very acute and precise about the results of the sum over the polar- 
ization I Y^ s . \Mjf\ 2 , we devote a whole section to the calculations of the various traces 
that intervene in the formal expression of the unpolarized DCS given by Eq. (0). We have 
to calculate : 

^EEK^i 2 = ^EEi^/' s /)[^ + c ^ + ^2^fe, Sl )i 2 

Si Sf Si Sf 

= 5££|S(P/>*/)A ( *Vpi,«0f (35) 

Si Sf 
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with : 

A« = [ 7 ° - 2k a 2 MPiHPf)]B s 

+ [c( Pi )-f ^ 2 + c(p f )4 2 ^ ]B 2s (36) 
using standard techniques of the 7 matrix algebra, one has : 

^EEK?! 2 = ^r{(^c + c 2 )AW(^c + c 2 )A (s) } (37) 

Si Sf 

with : 

a (s) = 7 o A wy 

= [7 - 2k a 2 MPi)c(PfW s 
+ [cfe)^ 1 ^7 + c(p / )7 #i]5L 

+ HpMl + c( Pf h°mB* 2s (38) 
There are nine main traces to be calculated. We write them explicitly : 

M l = Tr{{i> f c + c 2 )C (frc + c 2 )C }\B s \ 2 
M 2 = Tr{{i> f c + c 2 )C (Ac + c 2 )^ 
M 3 = Tr{{j> f c + c 2 )C (Ac + c 2 )C 2 }B s B* s 
M A = Tr{{j)fc + c 2 )d(Ac + c 2 )C }B* s B ls 

M 5 = Tr{(j> f c + c 2 )C 1 (frc + c 2 )C 1 }\B ls \ 2 (39) 
M 6 = Tr{{j) } c + c 2 )d(Ac + c 2 )C 2 }S 1SJ B 2 * s 
M 7 = Tr{(^ /C + c 2 )C 2 (Ac + c 2 )C }B 2s B; 
Ms = Tr{{i> f c + c 2 )C 2 (frc + c 2 )^}^^ 
M 9 = Tr{^ f c + c 2 )C 2 ^ t c + c 2 )C 2 }\B 2s \ 2 

To simplify the notations, we will drop the argument of the various ordinary Bessel functions 
that appear. The diagonal terms give rise to : 

Mi oc \B S \ 2 = J 2 

M 5 oc |5 ls | 2 = \{ J 2 +1 + 2 J s+1 J s _! cos(20 o ) + J 2 _J (40) 
M 9 oc |5 2s | 2 = i(J s 2 +1 - 2J S+1 J s _! cos(20 o ) + J 2 ^) 

So, taking into account the fact that the traces multiplying \B S \ 2 , \B ls \ 2 and |-B 2s | 2 are not 
zero, one expects that terms proportional to J s+ i J s _i cos(20 o ) will be present in the expression 
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of the scattering DCS. The first controversy between our work and the result of Szymanowski 
et al M concerns the traces Aia ancLMs- Since : 



M G oc B ls B* s = f(J 2 +1 - 2iJ s+1 J s _! sin(20 o ) - J 2 ^) 
M 8 oc B* u B 2s = f ( J s 2 +1 + 2iJ s+1 J s „ 1 sin(20 o ) - J a 2 _i) 



(41) 



and with little familiarity with the 7 matrix algebra, one can see at once that if the cor- 
responding traces are not zero then the net contribution of M.§ + Ais will contain a term 
proportional to J s+ i J s ~\ sin(20o)- We shall demonstrate that in what follows. We have : 

M s = Tr{{p} f c + c 2 )C 1 (frc + c 2 )C 2 }B ls B* 2s 

= Tr{(f> f c + c 2 ) [c( Pl h°Ui + c(p f )fah°KtiC + c 2 ) 

[c(p<)^7° + cip^U^BuB^ (42) 

From now on, we define a 4-vector : 

7f = (1,0,0,0) (43) 



We can therefore write : 



Then, Eq. (41) becomes : 



l=il (44) 



M 6 = Tr{(j> f c+ c 2 )C 1 (frc + c 2 )C 2 }B ls B* 2s 

= TrUp-fC + c 2 ) [cfaMfa + c(p f )fiiM(& c + c 2 ) 

[cipMil + cipfWfcftBuBl (45) 

In M, the authors claim that the controversial sin(20 o ) term disappear because it is propor- 
tional to terms like Tr{(^/c + c 2 )7°^i(^c + c 2 )^2^7°}- This term as well as Tr{{j[>fC + 
c 2 )^i^7°(^jC + c 2 )7°^2} are indeed zero but for Tr{{j)fC + c 2 )7°^i(^.jC + c 2 )7°^2} and 
Tr{(J)fC + c 2 )^i^7°(^jC + c 2 )^2^7°} this is no longer true. These terms are not zero and 
we give explicitly their values : 

Tr{(^c + c 2 )7 ^i(^c + c 2 ) 7 ^2} = Tr{(^ f c + c 2 )^ 7 (/^ + c 2 )^7°} 

= %w 2 {(a 1 .p f )(a 2 .p i ) + (a 1 .p i )(a 2 .pf)} (46) 

In most case, the various traces are zero except when the cyclic process of taking scalar 
products of pairs comes to products such that : 

(k.TjXk.rjXa^PiXa^pf) _ 
(A;.7/)(A;.^)(ai.p/)(a2.pi) 
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in which case, one has contributions proportional to w 2 (a>i.pi)(a 2 .pf) and w 2 (ai.pf)(a 2 .pi) 
respectively. Explicitly, we give the result for A4q and Ais- One has : 



w 

— {2sin(20 o ) 

cr 



(ai-Pi) (ga-P/) + (a 2 -P») (ai-Pf) 
_ (k.pi) (k.pf) (k.pi) (k.p f ) 



Js+lJs-l 



■{(a 1 .p i )(a 2 .pf) + {ai.pf){a 2 .p i )}J 2 s _ 1 
+{(a 1 .p i )(a 2 .p f ) + (a l .p f )(a 2 .p i )}J 2 +l }} 



(48) 



while A4s is given by : 



w 



{2 sin(20 c 



(ai-Pi) (og-P/) (a 2 -Pi) (ai-P/) ' 
_ (A;.;?*) (fc.p/) {k.pi) (k.p f ) 



Js+lJs-l 



-i[-{(a 1 .p i )(a 2 .Pf) + (ai.p / )(a 2 .pi)}J s 2 _ 1 
+{(ai.pi)(o2.p/) + (ai.p / )(a 2 .p i )}J s 2 +1 ]} 



(49) 



The fact that complex numbers appear in the expressions of A4q and A4s is not surprising 
since the former is the complex conjugate of the latter and their real sum is such that : 



4w^ 



sin(20 o ) 



(ai-Pi) (a 2 .p f ) (a 2 .pi) {ai.pf) 



Js+lJs-l 



(50) 



(k.pi) (k.p f ) (k.pi) (k.pf) 

So, the first controversy is settled and there is indeed a term containing sin(20o) in the 
expression of the scattering cross section. To put an end to any further criticism, we give 
in the Appendix the Reduce program we have written with the necessary commentaries and 
observations so that anyone in the scientific community having some knowledge of this powerful 
symbolic computational software can easily try it and check our results. Before writing our 
Reduce program, we have extensively studied the textbook by A. G. Grozin || which is full 
of worked examples in various fields of physics particularly in QED. We give the final result 
for the unpolarized DCS for the Mott scattering of a Dirac-Volkov electron : 



Z 2 19/ 



■ X 



jz{J 2 s A + (J s 2 +1 + J^)B + (i+*J*l x )C+ J S (J S _! + J S+1 )D} 



(51) 



where for notational simplicity we have dropped the argument z in the various ordinary Bessel 
functions. The coefficients A, B, C and D are respectively given by : 



A = c i -(q f .q i )c 2 + 2Q f Q i - 



(fc-g/) j_ ( k -<li 
(k.qi] 



a 2 u 2 



(a 2 ) 2 uj 2 



c 4 (k.q f )(k.qi 



+ 



a 2 uj 



(Qf - Qi) 



(k.q f )J c 2 (k.q f )(k.qi 
1 



■((Qf-Qi 



c 2 ) + 



(k.qi) (k.q f ) 



(52) 
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(a 2 ) 2 u; 2 



B 



+ 



OJ 



(ai.g/) (oi.gi) (o 2 .g/) («2-% 



2c*(k.q f )(k. qi ) 2c 2 V (*■?/) (fc-ft) (*■?/) (Mi) 



a 2 



a 2 ( (k.q f ) (k.qj) . 
4 l (^) (M/) ■ 



a 2 u; 2 



2c 2 



2c 2 (k.q f )(k.qi 
1 



■((?/•*) -c 2 ) + 



(fc.g/) (fe.gi 



(53) 



C 



a; 



\k.q f )(k.qi 



(cos(20 o ){(ai.g/)(a 1 .gi) - (a 2 .g/)(a2-9»)} + 



sin(20 o ){(ai.g/)(a 2 .ft) + (ai.Oj)(a 2 . <?/)}) 



(54) 



* - |«^) + (^))-|(g(A ft ) + g(A« f )) + 

w fQi(A.g f ) | Q/(A. 9i ) \ (55) 



c V (fc.g/) (fc.g- 
where A = a\ cos(0 o ) + a 2 sin(0 o ). 



3.1 Comparison of the coefficients. 

The argument about the missing term proportional to sin(0 o ) having been given a convincing 
explanation, we now turn to other remarks along the same lines since there are indeed other 
differences between our result and the result of [|l[] • We discuss now the difference occurring in 
our expression of the coefficient A and the corresponding one of [|lj. To make the comparison 
easier we give explicitly the simple relations between our coefficients and the corresponding 
coefficients of 0. One has : A(0J) = A/c 2 , 5(0) = 2£/c 2 , C(§) = 2C/c 2 and L>([JTJ) = 
2D/c 2 . In their expression multiplying the product 2J 2 (£), the single term ^ should 

come with a coefficient ~. In the appendix, we give a second Reduce program that allows the 
comparison between the coefficient A of fly and the coefficient A of this work. There are so 
many differences between our result and the result they found for the coefficient B that we refer 
the reader to our main Reduce program. The coefficient C has already been discussed. As for 
the coefficient D, we have found an expression that is linear in the electromagnetic potential. 
In the appendix, we give a third Reduce program. It is shown explicitly that if we ignore the 
first term in the coefficient multiplying J s (J s _i + J s +i) given in [[J, one easily gets the result 
we have obtained. This term does not come from the passage from the variables (p, p) to the 
variable (q, q) . The introduction of such 4- vector q is not useful, makes the calculations rather 
lengthy and gives rise to complicated expressions. As a supplementary consistency check of 
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our procedure used in writing the main Reduce program, we have reproduced the result of 
the DCS corresponding to the Compton scattering in an intense electromagnetic field given 
by Berestetzkii, Lifshitz and Pitaaevskii M. 

4 Conclusion. 

In this comment, we derived the correct expression of the first Born differential cross section 
for the scattering of the Dirac-Volkov electron by a Coulomb potential of a nucleus in the 
presence of a strong laser field. We have given the correct relativistic generalization of the 
Bunkin and Fedorov treatment that is valid for an arbitrary geometry. To prove that our 
results are correct, we give the Reduce program to let the scientific community judge their 
accuracy. We are adamant that the core of the whole controversy stems from the fact that 
in [|TJ, the vector introduced in Eq. ( fffi) of our work has not been properly dealt with 
while it is the common method to use when a trace contains a 7 matrix. Any standard QED 
textbook introduces this very elementary method. 

Appendix 

We give the main Reduce program that calculates the traces in Eq. fl3?|) . For this program to 
be readable, before every line, we give an explanation of the different instructions. Some are 
obvious, other are less straightforward and we also give the number of the equation to which 
it refers in the text wherever it is possible. In a Reduce program, a commentary is preceded 
by the symbol %. A Reduce instruction is not preceded by any symbol. 

The main program 

% This program calculates the trace appearing in Eq. (p7|) of the text. 

% The result must be multiplied by 4. Reduce calculates the quarter of any trace. 

% This is well explained in the manual ||]. 

% We first define the vector pf, pi, a, ai, a^, k, 77 (stands for nu), q^ qf. 
vector pfin, pin, aa, al, a2, k, nu, qin, qfin; 

% The command mass associates the relevant scalar variable as a mass with. 

% The corresponding vector. In the next instruction, cv stands for the velocity of light. 

mass k=0, pfin=cv, pin=cv; 

% The command mshell put a particle 'on the mass shell'. 
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% A substitution <vector variable>.<vector variable>=<mass>**2 is set up. 
mshell k, pun, pin; 

% The results of the above instruction are : k 2 = 0, p 2 = c 2 and pf = c 2 . 
on div; 

% We define properties of the vector 77 introduced in Eq. (|43|). 
let nu.nu=l, nu.k=w/cv, k.nu=w/cv; 

let nu.pfin=efin/cv, pfin.nu=efin/cv, nu.pin=ein/cv, pin.nu=ein/cv; 
let nu.al=0, al.nu=0, nu.a2=0, a2.nu=0; 
% We define the Maxwell gauge condition, 
let al.k=0, k.al=0, a2.k=0, k.a2=0; 

% We define the properties of the electromagnetic field potential. 

% We cannot use the variable a in our program. 

% Reduce interprets it as the matrix 7 5 . 

let al.al=aa.aa, a2.a2=aa.aa, al.a2=0, a2.al=0; 

% We define (pc + c 2 ). The variable '1' denote the fermionic line. 

for all p let gp(p)=cv*g(l,p)+cv**2; 

% We define the properties of the various quantities stemming from Eq. (|24]). 

% impart denotes the imaginary part and repart denotes the real part. 

% js, jspl and jsml denote respectively J s , J s+ i and J 8 -i- 

let impart(js)=0, impart(jspl)=0, impart(jsml)=0; 

let repart (js)=js, repart (jspl)=jspl, repart (jsml)=jsml; 

let impart(s)=0, impart(phiO)=0; 

let repart (s)=s, repart (phi0)=phi0; 

% We define the quantities of Eq. (p^). 

bs:=js*exp(i*s*phiO); 

bls: = (jspl*exp(i*(s+l)*phi0)+jsml*exp(i*(s-l)*phi0))/2; 
b2s: = (jspl*exp(i*(s+l)*phi0)-jsml*exp(i*(s-l)*phi0))/(2*i); 
repart (bs):=js*cos(s*phiO); 
impart (bs) : =j s*sin(s*phiO) ; 

repart (bis) : = (jspl*cos((s+l)*phi0)+jsml*cos((s-l)*phi0))/2; 

impart (bis) :=(jspl*sin((s+l)*phi0)+jsml*sin((s-l)*phi0))/2; 

repart (b2s): = (jspl*sin((s+l)*phi0)-jsml*sin((s-l)*phi0))/2; 

impart (b2s):=(jsml*cos((s-l)*phi0)-jspl*cos((s+l)*phi0))/2; 

% We ask Reduce not to perform the various traces for the time being. 



13 



nospur 1; 

% We define the various products appearing in Eq. (|37|). 
tl:=gp(pfin); 

t2:=(g(l,nu)-2*c*cp*(aa.aa)*(k.nu)*g(l,k))*bs; 

t3: = (c*g(l,nu)*g(l,k)*g(l,al)+cp*g(l,al)*g(l,k)*g(l,nu))*bls; 

t4:=(c*g(l,nu)*g(l J k)*g(l,a2)+cp*g(l,a2)*g(l,k)*g(l,nu))*b2s; 

t5:=gp(pin); 

t6: = (g(l,nu)-2*c*cp*(aa.aa)*(k.nu)*g(l,k))*conj(bs); 
t7:=(cp*g(l,nu)*g(l,k)*g(l,al)+c*g(l,al)*g(l,k)*g(l,nu))*conj(bls); 
t8: = (cp*g(l,nu)*g(l,k)*g(l,a2)+c*g(l,a2)*g(l,k)*g(l,nu))*conj(b2s); 
% To obtain compact expressions, we define : 
for all fi, let cos(fi)+i*sin(fi)=exp(i*fi); 
for all fi, let cos(fi)-i*sin(fi)exp(-i*fi); 

% We explicitly define the coeficients c(pi) = c{qi) and c(p/) = c(qf). 

c:=l/(2*cv*(k.qin); 

cp:=l/(2*cv*(k.qfin); 

% We define the product J s (J s _i + J s +i) of Bessel functions, 
let js*jsml+js*jspl=jsom; 

% We define the relation between the 4-vector p and q. 

let pin=qin+(aa.aa)*c*k/cv; 

let pfin=qfin+(aa.aa)*cp*k/cv; 

% Same definition for the energy. 

let ein=gqi+(aa.aa)*c*w/cv; 

let efin=gqf+(aa.aa)*cp*w/cv; 

% We load the package ASSIST [|J that simplifies the calculations. 
load_package assist; 

% We now ask Reduce to calculate the traces related to the fermionic line T. 
spur 1; 

% Each trace resz is calculated separately corresponds to the trace M.i of Eq. (|39|). 

resl:=tl*t2*t5*t6; 

trigreduce resl; 

res2:=tl*t2*t5*t7; 

trigreduce res2; 

res3:=tl*t2*t5*t8; 
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trigreduce res3; 
res4:=tl*t3*t5*t6; 
trigreduce res4; 
res5:=tl*t3*t5*t7; 
trigreduce res5; 
res6:=tl*t3*t5*t8; 
trigreduce res6; 
res7:=tl*t4*t5*t6; 
trigreduce res7; 
res8:=tl*t4*t5*t7; 
trigreduce res8; 
res9:=tl*t4*t5*t8; 
trigreduce res9; 

% The total trace is the sum of all these traces. 
restot:=resl+res2+res3+res4+res5+res6+res7+res8+res9; 
trigreduce restot; 

% The program is complete and Reduce takes seconds to give the answer. 

The most controversial terms correspond to res6 and res8. It is very easy to check that our 
claim is well founded and that there is indeed a term proportional to sm(20 o ) that is missing 
in We now turn to the Reduce program that shows that the single term ^ ^ should 
come with a coefficient | 

The coefficient A. 

This coefficient comes from the calculation of M.\ in Eq. (|39|). The coefficient A of IQ. is 
transformed so that no vector q appears anymore. The symbols qit and qft stand for q~i and 
q~f respectively. The instruction fl:=cl is the correct instruction that allows to find exactly 
the same result as that given in Eq. (|52"|). 

on div; 

vector qi, qf, qit, qft, k; 
al:=c**2-qi.qf+2*gqi*gqf/c**2; 
bl:=aa**2*(k.qft/k.qi+k.qit/k.qf)/(2*c**2); 
cl:=aa**4*w**2/(2*c**6*(k.qi)*(k.qf)); 
dl: = l-aa**2*w**2/(c**4*(k.qi)*(k.qf)); 
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el:=-aa**2*(l-k.qft*k.qit/(k.qi*k.qf))/(2*c**2); 

% Here, we use the fact that there is a factor 1/2 in the above mentioned single term. 
fl:=cl; 

gl:=aa**4*w**2*(k.qf*k.qft+k.qi*k.qit+aa**2*w**2/c**4); 

hl:=gl/(2*c**6*(k.qi)**2*(k.qf)**2); 

k.qit:=-k.qi+2*w*gqi/c**2; 

k.qft:=-k.qf+2*w*gqf/c**2; 

res:=(al+bl+cl)*dl+el+fl+hl; 

This program gives exactly our coefficient A given in Eq. ( |52"| ) 

The coefficient D. 

For the coefficient D of , we give the Reduce program omitting the first term that contains 
a factor that is quadratic in a. Doing so, we find exactly the same expression as that given in 
Eq. ©. 
on div; 

vector qi, qf, qit, qft, k, aron; 
terml:=0; 

term2:=aron.qi*k.qit/(c*k.qf)+aron.qi*k.qft/(c*k.qi); 

term3:=aron.qi/c+aron.qf/c; 

k.qit:=-k.qi+2*w*gqi/c**2; 

k.qft:=-k.qf+2*w*gqf/c**2; 

resd:=terml+term2+term3; 

We have given convincing arguments to support our results. We also gave programs that 
will allow anyone to check every stage of our reasoning and to reach the same conclusion : 
there are indeed mistakes in [|IJ. 
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